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The research discussed in this report was accomplished as part of 
the Maintenance Effectiveness Engineering Graduate Program conducted 
jointly by the DARCOM Intern Training Center and Texas A&M University. 

“AS such; the jdeas, concepts and results herein presented are those of 
the author and do not necessarily reflect approval or acceptance by the 
Department of the Arny. 
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ABSTRACT 


Research Performed bv Craig D. Hunter 
Under the Supervision of Dr. S. Bart Chiids 


Regression analysis normally implies the use of algebraic equations 
to describe a system; however, some cases would better be modeled by 
differential equations. This is accomplished by assuming a differential 
equation model for a given set of data and estimating the values of the 
unknown parameters within the model. These values are then systematically 
perturbed to generate particular solutions which are superimposed to yield 
a better estimate of the unknowns. This process is repeated until a 


specified accuracy is met. 
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Through an analysis of variance, the statistical characteristics of 


linear regression can be generated for most nth 


order differential equa- 
tions. This provides a basis for evaluating the ‘acceptance or rejection’ ; 
of the regression. 

The characteristics generated consist of an ANOVA table (uncorrected), 
general F test on the regression, the pe value, covariance matrix of the 
superposition constants, an estimate of the variance about the regression, 


an estimate of the variance of the parameters, and the confidence inter- 


vals on these estimates. _ 
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Through an analysis of variance, the statistical characteristics of linear 
regression can be generated for most nth order differential equations. This 
provides a basis for evaluating the ‘acceptance or rejection! of the reqression 


The characteristics generated consist of an ANOVA table (uncorrected), 
general F test on the regression, the R* value, covariance matrix of the super- 
position constants, an estimate of the variance about the regression, an esti- 


mate of the variance of the parameters, and the confidence intervals on these 
me estimates, 
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CHAPTER I 


INTRODUCTION 


Regression analysis usually impiies determining a set of independent 


variables and corresponding coefficients for algebraic equations such 
that these equations best portray the data or system given. Examples 
are presented in which algebraic models are inadequate and can be appro- 
priately modeled by differential equations. 

This report is an investigation into the validity and accuracy of 
a differential equation regression package, both linear and non-linear 
(Childs, 1970). It is not the intent to deal with extensions of the 
established theory behind such an approach (Bellman and Kalaba, 1965), 
(Roberts and Shipman, 1972). Nomenclature and background material are 
given in Chapter II. 

A control program was written to provide a means for evaluating 
QUASII (Childs, 1970). The control program yields solutions of a speci- 
fied accuracy (within the computational restrictions of the computer) 
providing a means for comparison. QUASII would be used instead of a 
program yielding "exact" solutions because of ease of use, 

The control program is based on power series evaluation and inte- 
gration of differential equations; whereas, QUASII utilizes a 4th order 
Runge-Kutta integrator. The Runge-Kutta integrator allows simpler pro- 
gramming of non-linear differential equations, 

Two solution areas are of interest: 

Case I: The number of boundary conditions (m) equals 
the order of the diff2rential equation (n), 


therefore all boundary cc.iditions are to be 


met exactly (exact system). 

Case II: The number of boundary conditions (m) exceeds 
the order of the differential equation (n), 
therefore some boundary conditions are to be met 
exactly and others to be met in a least squares 


sense (over-determined system). 
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Both QUASII and the control program utilize a least squares fit of 
the non-exact boundary conditions in Case II. Further study in this area 


is being done by Walker (Ref. 21). These results will be available in 
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April 1976 and entails meeting these boundary conditions in a least 
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squares and other criteria. 


Case I is of theoretical interest. This report contains results for 
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Case II studies. 
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CHAPTER II 
DIFFERENTIAL REGRESSION 


The following example of a linear case from Childs et. al. 1971, 
(Ref. 6), is for clarification of terms and background. For information 
on non-linear cases, refer to Reference 5 and 17. 

Given: A simple mass-spring-damper system governed by the equation; 

x + ux + Ex = sin(wt) (2-1) 
where (+) denotes the total derivative of that term with respect to t, 
the independent variable. 


Denoting the following state variables; 


yy = xX Yo = xX (2-2) 
equation 2-1 can be rewritten as two non-trivial first order differential 
equations: 

41° V2 
; (2-3) 
Yo = -EY) ~v¥y + sin(wt) 
The state variables y, and Yo are functions of t, y,(t) (t has been omit- 
ted for simplicity). 


If wand — are known, equations 2-1 and 2-3 are linear differential 
equations. Assuming that p and — are unknown, which is not an unreason- 
able assumption, two more state variables must be added; 

¥3 = Yrs 
ot which the fcllowing four equations result; 


yz = 0 (2-4) 


ti 


Y4 0 (2-5) 
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yy = Yo (2-6) 


Yo = “Yy¥q ~Yo3 * sin(wt) (2-7) 


the frequency w could have been assumed unknown, introducing a fifth 
state variable, Ye = We Except for equation 2-7, the four differential 
equations are linear eith 2-4 and 2-5 being trivial. Since equation 2-7 
is non-linear, the system is non-linear. 

At least two Newton type or Taylor series methods used in solving 
such systems are perturbation and quasilinearization methods. 

Perturbation methods are used in the control program and QUASII for 
ease of programming. For information regarding this theory, see Refer- 
ence 5, Doiron (1970), or any of Childs' writings cn non-linear differen- 
tial equation solutions. 

The right hana side (RHS) of equation 2-3 is unknown (except for yu 
and &), therefore, avenues are opened for the implementation of boundary 
conditions or constraints which must be met. In a mass-spring-damper 
system governed by equation 2-1, the use of an accelerometer to measure 
the acceleration of the mass at certain points in time (t,) is within 
experimental procedure. An accumulation of constraints (boundary values) 
could be obtained in the following form; 


x(t.) = b. i= 1...m (2-8) 


where b. is the measured acceleration at tine t, and m is the number 
of boundary values. The LHS of equation 2-3, Yo (yp = x), can now be re- 
placed by bis 
c ‘ _q\ 
bs = -€y,(t,) -uyp(t.) + sin(wt,) (2-9) 


in more general terms; 


q, (y(t;)) # b, i= l...m (2-10) 


where CF is a boundary condition operator. Equation 2-10 defines m 
boundary conditions resulting in m linear (or as in equation 2-7, non- 
linear) equations. From Chapter I, m must be equal to or greater than n, 
the order of the differential equation (number of state variables). These 
m differential equations will either be met exactly or at least in a least 
square (best fit) manner. 


A brief summary of a "shooting" technique of solving linear differen- 
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tial equations follows. Given; 
y=lyt+f (2-11) | 
where L is a linear operator and f is a forcing function. 

A set of homogeneous solutions is usually superimposed upon a par- 
ticular solution. The number of homogeneous solutions is equal to n, the 
order of the differential equation (see equation 2-1, n=2). Assuming 
that no initial conaitions are known, n = r (let r equal the number of 
unknown initial conditions, where r <n). When only r initial conditions 
are unknown, r homogeneous equations are required. The result of the 


superposition of a particular solution (2-12) and a set of homogeneous 


equations (2-13) is: 


p=lp+f (2-12) 

A= LH (2-13) 
t p(k) 

y=pt+He=p+ > h Bk : (2-14) 


It is an established theory that the sum of a particular solution 
and a set cf homogeneous equations yields another particular solution, 


An alternative is to superimpose only particular solutions. This scheme 


is used in QUASII and the control program; 


y= Pas p(k) a, (2-15) 


where 

p(k) zs Lp(k) +f (2-16) 
There will be r+] particular solutions, where r is the number of unknown 
initial values, 


A restriction on the superposition constants is: 
r 
L ay = 1 . (2-17) 


This is a result of equation 2-11, where the differential equation is 
given by: 
yely+f 
If particular solutions of the form of equation 2-16 are superimposed 
without the constraint 2-17, the resulting differential equation could be 
y= Ly + Cf (2-18) 


where 
} (2-19) 
c= a, #1 2-19 
K=o * 


The original differential equation has not been satisfied and necessi- 


tates the restriction on the superposition constants to yield 


r 
ye lyst) ap | (2-20) 


The computational strategy requires choosing initial arbitrary values 
for the r unknown state variables. These values are used in defining the 
unperturbed particular solution (denoted by p')), Each of the unknown 


variables (r of them) will in turn be individually perturbed by a 


os 
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predetermined (aribtrary) constant to yield an appropriately perturbed 
particular solution. These particular solutions, obtained by integration, 
yield values corresponding to the boundary values at the appropriate value 
of the independent variable. 

The result is a matrix equation of the form 
Sa=d (2-21) 
where S is a function of the integrated solutions, a is the vector of 
superposition constants, and d is a vector of boundary values. The re- 
sulting m+1 equations (m due to the boundary conditions and one due to the 
restriction on the superposition constants) contain r+1 unknown a's 
(r+1 < m+1), which can be determined as discussed in Chapter III. 

The resulting y solutions (state variables) are the initial arbi- 
trary estimates plus the addition of the corresponding superposition 
constant times the perturbation constant. 

Evaluation of the a's was stated to be a rather simple process. This 
entails a Gauss-Jordan reduction utilizing a least squares fit of the 
over-determined system. A more in-denth discussion of the theory is pre- 


sented in Chapter III and Reference 4. 
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CHAPTER III 
THE OVER-DETERMINED SYSTEM 


When the number of boundary conditions (m) equals the order (n) 
of the differential equation, the boundary conditions can theoretically 
all be met exactly, however, in the case of the over-determined system, 
this is impossible. Thus, several boundary values will have to be met 
by some other means. 

The problem of Chapter II was narrowed to the final equation: 

Sa=d (3-1) 

This equation and the role.of a least squares fit are the topics of this 
chapter. 


Examining the composition of the S matrix, two subdivisions can be 


shown; 
Si k+l 
Sp m-k 
rt+i 


where k is the number of boundary conditions to be met exactly. The Sa 
matrix contains the exact conditions and SB the least squares boundary 
conditions, | 

A further reduction or partitioning can be; 


Saye we al = |d (3-2) 


ommngerene 


where 
S) is order (k+1) x (k+1) 
So is order (k+1) x (r-k) 
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S, is order (m-k) x (k+1) 
Sy is order (m-k) x (r-k) 


By a Gauss-Jordan reduction, the Sy matrix can be transformed into 


an identity matrix resulting in S3 being a null matrix. The results are; 


t 
I $5 an de 
mee ees | eee (3-3) 


where the subscripts 1 and e refer to least square and exact boundary 
conditions respectfully. 
The null matrix (S3) allows for an equation (3-4) containing only 
least squares conditons. 
S44 = dy (3-4) 
Equation 3-4 yields m-k equations each with an unknown. 
Premultiplying by the transpose of S4 gives the normal equation 


(refer to Appendix B); 

(S4)"sta, = (Sp)'as (3-5) 
resulting in the solution for ay: 

ay = [(spr'stg} sy)" (3-6) 
Having solved for a, in equation 3-6, a, (from exact conditions) can be 
solved by referring to equation 3-3; 

la, + $58) = de (3-7) 
or 

= qd! i" 

Equation 3-4 and 3-5 are important and will be used extensively in 


Chapter IV. 
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CHAPTER IV 


ANALYSIS OF VARIANCE 


The boundary conditions are: 


a(y(t.)) = b, 12 decom (4-1) 
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This defines an operator q operating on the y(t.)'s to yield the desired 
boundary values. The statistical basis for ‘accepting or rejecting’ the 
regression is the nearness of the equality 4-1. 
In Appendix B, the equation; 

z=Wdte (4-2) 
is analyzed through an analysis of variance. 


From equations 4-1 and 4-2, the following results can be obtained; 


b=qte (4-3) 
where 

q = aly) 

b=2 

€=e 

q = Nd 


Note that q(y(t.)) may not be linear as required for an analysis of 
variance as described in Appendix B, but the solution process is itera- 
tive using linearized equations, Because the steps aire small, the assump- 
tion of q being linear is a good approximation. 


An analysis of variance can be accomplished as shown in Table 4.1. 


as LEE Ee OS Seek ce eretag 
sais essai tS BRE se ta ae NL MRE “aati ts 
; : 


mete Te ae 


PRR re vettattonnnrrn ae» 


| | n 
if TABLE 4.1 
cas 
i} 
: i ANOVA TABLE 
? Sum of Degrees of 
f Source Squares Freedom Mean Square 
an 
: | Due to regression q'b r SS/r 
a: About the regression b'b - q'b m-k-r 5 
3 (residual) 
4 
f Total (uncorrected) b'b m-k 
‘ 4 
‘ The following are calculated directly from Table 4.1: 
i 
me se SSregression - (m-k)(mean of by's)? (4-4) 
SStotal - (m-k)(mean of b,'s) 2 
ss - MS_regression (4-5) 
MS residual 
sé = MS residual = estimated variance of system (4-6) 


The Rong value will be compared to a Fisher's F with: 


Probability of lea (a is the producer's risk) 
Numerator degrees of freedom r 
Denominator degrees of freedom (m-k-r) 

If F eal exceeds Table F, the regression is 'accepted'. (Note: The 
producer's risk is defined as the probability of rejecting the null 
hypothesis (Hy) (see Appendix B) when in fact the null hypothesis is true. 

Confidence limits for the estimated least square boundary values can 
be determined once the individual variances are known. Referring to 
Appendix B, the variance of the estimated observation is given by equa- 


tion B-10: 


: 
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est. var(z,) = Wl (WW) "Ws (4-7) 


Equation 4-3 shows that q is synonymous with Wd of Appendix B, but q 
cannot be subdivided into a matrix of independent variables and a vector 
of parameters, thcrefore another approach must be used. 

Using the system of equation 3-5, a relation between d, (least 
squares boundary values of b) and dy (transformation ot dy) can be deter- 
mined, 


The variances of the individual least squares boundary values; 


él] Lf Ms 


is equal to the variances of the transformed system: 
i) 
o si [+] = [4 (4-9) 
: L 


var (d, ) = var (d> ) (4-10) 
The variances of the estimated least squares boundary values are 
given by; 
est. var(b,) = Sy, Lisg)'sg] Psy y's"? i=1...m (4-11) 
where 54, is the row vector of S4 corresponding to the jth least squars 


boundary condition, 


The resulting confidence limits are; 


b, # t(v,l- 3) [est. var(6,)] (4-12) 
where 

v = m-k-r (4-13) 
i= k+l ... m 


A main objective of the regression package is to obtain the solutions 
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for the y;'s given in Chapter II. Table 4.1 does not always provide a 
means for irvestigzting this because it involves the boundary condition 
functions. By doing a partial analysis of variance on equation 4-10, 


the variar.ce of the a's can be determined; 
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est. var(a,) = ((st)isty7tg!2 (4-14) 
7] 4’ ~4 
: where 
rf ry Teas - Trerylan 
| Pee ait al (4-15) 
m-k-r 


Referring to Chapter II, the ys solutions are equal to the previous 
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value (initial guess) plus the perturbation times a superposition con- 


stant; 


a a 


Y * Yorevious + (a, * perturbation) (4-16) 
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therefore 


a ae 


var(y) = var(a,) * (perturbation)* (4-17) 


This results in a confidence interval for these solutions as follows: 


a aed 


y,tt(vyl- 3) [est. var(9,)| oe aes y (4-18) 
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CHAPTER V 


RESULTS-COASTING DATA 


The statistical study outlined in Chapter IV was accomplished by a 


subroutine package PSTAT (Appendix D contains PSTAT and supportive sub- a4 

routines). ae 
Since experimental data was not available on the system of equation 

2-1, a problem of determining the aerodynamic drag coefficient and roll- 

ing friction coefficient of an automobile was used (Ref. 21). The non- 

linear equation governing this system is; 


ee A p 02 
x + OFF OC, x + weg = 0 (5-1) 
i 


where 


Be A ane He 


o = air density (slugs/ft?) 
A= frontal area of vehicle (Ft?) 


Tee Be 


M= mass of vehicle (slugs) 
g = acceleration of gravity (ft/sec”) 
Cy coefficient of drag 
ue rolling friction coefficient 
x = velocity (ft/sec) 
x = acceleration (ft/sec?) 

This problem demonstrates the advantages of modeling by differential 
equations with a regression package like QUASII. Attempts to determine 
Cy and We by other than numerical methods (such as wind tunnel testing, 
treadmills, etc.) would be more laborious, less accurate, and more ex- 
pensive. 


Coasting data (value of x) for a Sunbeam Alpine was obtained from 
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Road and Track Road Test Annual for 1966. Figure 5.1 is a listing of 
input for such a system. Figure 5.2 is an output of the statistical 
analysis from PSTAT. Figure 5.1 shows only one non-trivial differential 
equation. Although equation 5-1 is a 2nd order differential equation, 
it does not contain an x term, allowing a single first order representa- 
tion. 

Except for three parameters, the PSTAT output (Figure 5.2) has been 
previously defined. The three additional parameters are: 

1) Mean of the observations is the sum of the least squares observa- 
tions (boundary values) divided by the total number of these observations: 


m-k 
es ae (5.2) 
a a 

2) Coefficient of variation (cv) is a measure of the dispersion of 


the data (Hald, 1952) and is expressed as: 


_ MS residual 3 
~ mean of observations | s/b 


(5-3) 
The importance of this simple parameter is questionable. 

3) P(F(ALPHA).GY.FCAL) is the probability that a random variable 
that is F distributed, with the same degrees of freedom as F(ALPHA), will 
be greater than FCAL, thus the minimum producer's risk that can be assumed, 

The PSTAT output (Figure 5.2) is an ANOVA table of the uncorrected 
sum of squares plus the corresponding degrees of freedom and mean square 
terms, Next is R° and the coefficient of variation value. A Rr? value 
of 1 indicates a 'perfect' regression. The TEST OF OVERALL REGRESSION 
statement is a test of the hypothesis stated in Appendix B for a risk of 
ALPHA (a). If the FCAL exceeds F(ALPHA), H, is rejected and the alterna- 


tive, Hs is accepted. Under these conditions, 
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*eKKKACCEPT REGRESSION***** is printed. This is not to imply this is 


the 'best' possible regression, but for the stated hypothesis and the 
given data, it is an ‘acceptable’ regression. 

Following the hypothesis test are two sections showing confidence 
intervals on the observations and state variables respectfully. These 
limits are all based on the same producer's risk (ALPHA) and t value 
(Student's t distribution). In addition, the observed values (given 
boundary values) are shown with the estimated values from the regression. 
The residual terms; 

(observed value) - (estimated value) (5-4) 
and the sum of the differences are included. Both sections contain the 
estimated variance, or standard error, of the individually estimated 
solutions, 

The last section is the covariance matrix of the superposition con- 
stants (least squares). The diagonal elements are the variances of 
those constants, 

A regression cannot be 'accepted' or ‘rejected’ on the basis of 
just one statistic. The programmer must qualify an acceptable solution 
and evaluate the overall regression, confidence limits, variances, hypo- 
thesis, etc. for the given system to determine whether a resulting model 


is acceptable, PSTAT does furnish enough information to evaluate the 


regression validity in most cases. 
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CHAPTER VI 
CONCLUSIONS ON STATISTICAL RESULTS 


PSTAT, a subroutine from the control program, provides a means of 
evaluating the acceptability of differential equation regressions. The 
control program and QUASII are based on the same theory and procedures 
(except for the integration methods). PSTAT is presently not directly 
compatible with QUASII, but does provide a basis on which QUASII solutions 
can be evaluated. 

This statistical package can be used for linear and non-linear 
systems (see Reference 4 for types of non-linearities). 

PSTAT cannot handle multiple observations for distinct value(s) of 
the independent variable(s). The ANOVA table and procedures should be 
modified to include pure error and further hypothesis testing. 

The assumption of t'e linearity of equation 4-2 for the purpose of 
an analysis of variance, does not appear to be valid for examples worked 
with non-zero forcing functions (see equation 2-1). This area should be 
investigated from the possible approach uf a non-linear analysis of 
variance, 

In its present farm, PSTAT coupled with QUASII provides a powerful 
tool for statistically evaluating many systems not previously modeled 


in a sufficient manner by standard algebraic regression packages. 
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APPENDIX A 


NOMENCLATURE 


Frontal area of a vehicle. 

Superposition constant. 

Boundary condition. 

Coefficient of variation. 

Transformation of boundary condition. 

Expected value of, 

Forcing function, 

Theoretical value from Fisher's F Distribution Tables. 
Variable that is F distributed. 

Acceleration of gravity. 

Set of homogeneous solutions. 

Identity matrix. 

Number of boundary conditions to be met exactly. 
Linear operator. 

Mass. 

Number of boundary conditions. 


Mean-square. 


Order of the differential equation (number of state var.). 


A particular solution, 

Set of particular solutions. 

Operator, either linear or non-linear. 
Number of unknown initial conditions. 


Percent variation due to regression. 


Estimate of variance. 


Estimate of variance. 


Matrix of the integrated values of the particular solutions. 


_Sum of squares. 


Standard error of. 


Independent variable time or a student's t statistic. 


Variance of. 

State variable. 
Coefficient of damping. 
Spring constant, 
Frequency. 

Producer's risk, 
Superposition constant, 
Rolling friction coefficient. 
Air density. 
Coefficient of drag. 
Degrees of freedom. 
Estimated, 

Exact. 

Least squares. 


Transpose, 


Transformation of by matrix operations, 


Estimate. 


Mean. 
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APPENDIX B 


GENERAL LINEAR REGRESSION THEORY 
(Draper and Smith, 1966) 

To handle multiple linear regression, a matrix approach is desirable. 
Unless otherwise stated, all unsubscripted letters are matrices or vec- 
tors. 

Consider a linear problem (only in the coefficients) with n obser- 


vations and p independent variables, the following linear equation could 


arise; 
=Wet+e (B-1) 
where 
z is a (n x 1) vector of observations 
Wis a (n x p) matrix of independent variable 
8 is a (n x 1) vector of parameters to be estimated 
e is a (n x 1) vector of errors 


Certain basic assumptions must be made on the error terms, c; 


e,'S are random variables 


var(e;) = of (unknown) 


e, and ec, are uncorrelated (i#j) (B-2) 
cov(e; 5 .)=0 
e,'S are normally distributed; N(0,0°) 


e,'S are therefore independent 


As a result of cov(es,e5)=0, 2. and z, are also uncorrelated and 


By evaluating the error sum of squares and substituting the least 


squares estimates of g (d an n vector), the following normal equation 
results: 

(Wwf = we (B~3) 
Using the decomposition of the total sum of squares (uncorrected), the 


ANOVA (analysis of variance) Table B.1 results. 


TABLE B.1 

Source Sum of Squares Degrees of Freedom Mean Square 

Due to regression dWz p MSpeg : 
About regression zz - dy n=p MS = 5 ‘ 
(residual) i 

é 

Total (uncorrected) zz n 


MS, is an estimate of the variance about the regression. 

Once the ANOVA table is constructed, several conclusions can be 

deduced: 

I. A percent measurement of the total variation about the mean 
attributed to the regression (square of the multiple correlation 
coefficient) is given by; 

t, =e 


T 
Re S d W - NZ (B-4) 


zZ-Nnz 
where z is the mean cf the observations. 
A value of R? = 1 implies that the prediction is 'perfect'. 


II. To test the overall regression equation, the hypothesis: 


Hp! fy = 8, ee = 8, = 0 (B-5) 
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27 
must be tested against the alternate hypothesis: 
Hy: 8B. #0 1 = 0...p (B-6) 

This employs the Fisher's F distribution. The Foal statistic 
is: 

MS 

REG 
MSe 


This value is to be comnared to a table value of; 


Fe) 
where: 
l-a is a probability 
p is the degree of freedom for the numerator 
n-p is the degree of freedom for the denominator 
If Foal exceeds the table value, the null hypothesis is rejected 


implying that the variation in the data is more than would be 
expected by chance in (l-a)100% of similar data sets. 


The covariance matrix of the least squares estimates (d) is; 


est. var(d) = (Wu)74s2 (B-8) 
where the jth diagonal elements is the estimated variance of 
d.. Equation B-8 coupled with the student's t distribution 
yields a confidence interval for d. with a risk of alpha (a) of: 
d, + t(n-p, 1l-a/2) [estimate $.e@, (4,)| i=l...n  (B-9) 
Similarly, a confidence interval for the estimates of z (z) can 
be determined. Based on linear combinations of random variables, 


the estimated variance of z. is; 


est. var(z,) = W, (Wh) Wis? ieleean (B-10) 


where W, is the vector of W corresponding to the jth row or 


boundary condition (least square). The resulting confidence 


interval is: 


z. + t(n-p, l-a/2) fest. var(2,)| i=1...n (B-11) 
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APPENDIX C 


APPLICATIONS TO ALGEBRAIC REGRESSION 


Illustration of the versatility of QUASII and verification of the 
control program can be shown utilizing the following algebraic regres- 
sion: 


Bs 2 3 
x=batb, t+bot +bt (C-1) 


_ 


Equation C-1 is a linear third order regression equation, 


Utilizing the following state variables: 


3. 2 3 
yy RX=botd t+bot +bt 


yorieb,+2bpt+3bgt” (C~2) 


3 


Y3rx=2bo+6b4t 


dec ne RE Rr ont STE Se 


Yg=x=6b, 

Four differential equations result: 
Yi7V0 i 
YorY 
Y3Vq 
y4=0 

Equations C-3 are linear with Ya being trivial. 

By solving the set of differential equations C-3, the values of the 
y;'s at t=0 can be determined. Equation set C~2 yields the following 


results: 


b,=yp(0) ba=y4(0)/6 (C-4) 
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APPENDIX D 


THE CONTROL PROGRAM WITH OUTPUT OF THE 
NON-LINEAR MASS-SPRING-DAMPER SYSTEM 

The contre] program is based on the same theory as QUASII and both 
are quite similar in structure. Several subroutines are common (with 
slight differences) to both programs. However, QUASII is capable of 
handling a larger variety of problems. As stated earlier, the main 
difference is that QUASII utilizes a 4th order Runge-cutta integrator 
opposed to power series integration employed by the control program. Be- 
cause of this point a subroutine RECURL must be rewritten for each system 
to be modeled. 

The following pages contain a copy of the control program main-line, 
supporting subroutines, and post-statistical package for the non-linear 
system of equation 2-1. The output does not include the results from 
the PSTAT package. 


The program is well documented for both Case I and II problems, 
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APPENDIX D 


THE CONTROL PROGRAM WITH OUTPUT OF THE 
NON-LINEAR MASS-SPRING-DAMPER SYSTEM 

The control program is based on the same theory as QUASII and both 
are quite similar in structure. Several subroutines are common (with 
slight differences) to both programs. However, QUASII is capable of 
handling a larger variety of problems. As stated earlier, the main 
difference is that QUASII utilizes a ath order Runge-Kutta integrator 
opposed to power series integration employed by the control program. Be-~ 
cause of this point a subroutine RECURL must be rewritten for each system 
to be madeled. 

The following pages contain a copy of the control program main-line, 
supporting subroutines, and post-statistical package for the non-linear 
system of equation 2-1. The output does not include the results from 
the PSTAT package. 


The program is well documented for both Case I and 11 problems. 
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